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EXPERIMENTAL  INVESTIGATIONS  OW  THE  LIGHT  SCATTERING  OP  COLLOIDAL  SPHERES. 

V.  DETERMINATION  OP  SIZE  DISTRIBUTION  CURVES  BY  MEANS  OP  SPECTRA  OP  THE 
SCATTERING  RATIO1. 

Wilfried  Heller  end  M.  L.  Wailach2 
Department  of  Chemistry,  Wayne  State  University,  Detroit,  Michigan 

and 

A.  F.  Stevenson 

Department  of  Physics,  Wayne  State  University,  Detroit,  Michigan 

I.  INTRODUCTION 

The  preceding  paper  in  this  series  dealt  with  the  scattering 
ratio,  (J  ,  determined  at  an  angle  of  90°  with  respect  to  the  incident  beam 
as  a  means  for  determining  particle  site  from  light  scattering  in  mono- 
disperse  colloidal  c.isperaiona  of  non-absorbing  spheres  .  The  principal 
advantages  of  the  at  e  of  S1  over  a  series  of  other  possible  light  scatter¬ 
ing  techniques  were  found  to  be  fa)  its  relative  Insensitivity  to  errors 
in  concentration,  (V)  the  fact  that  no  instrument  constant  is  needed  and 
(c)  its  high  prospective  sensitivity  to  beterodlsperslon.  In  view  of  these 
advantages,  the  scattering  ratio  appeared  very  promising  as  an  argument  for 
determining  Rita  distribution  curves  in  bnterodisperse  systems.  An  exten- 

k 

slve  theoretical  investigation  of  this  problem  laid  the  foundation  for  a 
practical  application  of  <T -spectra  to  this  end.  The  present  paper  is 
concerned  with  this  practical  application^,  Its  results  apply  efually  to 
measurements  of  the  depolarization  (polarization  ratio  at  90°) • 


1.  This  work  was  supported  by  the  Office  of  Naval  Research.  The  results 
given  in  the  present  paper  were  presented  at  the  134th  meeting  of 
the  American  Chemical  Society;  Chicago,  September,  1958. 

2.  Present  addresa;  Film  Department.,  E.  I.  DuPont  &  Co.,  Wilmington, 
Delaware. 

3-  W,  Heller  and  R.  Tablblon,  J.  Phyn.  Chen.  66,  2039  (1962). 

4.  A.  F.  Stevenson,  W.  Heller,  and  M.  L.  Wailach,  J.  Cham*  Rqrs.  34, 

1789  (1961). 

3.  Regarding  other  recent  attempts  at  using  light  scattering ffor  a  deter¬ 
mination  of  size  distributions,  reference  may  be  made  to  the  extensive 
literature  review  given  in  reference 
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II.  BRIEJ'  REVIEW  01’  THE  THEORY 

g 

The  theoretical  vuriation  of  the  scattering  ratio  6 with  OC 

*7 

investigated  previously  shoves  (Figure  1  to  3  ia  ref*  7)  that  its 
spectrum  will  exhibit  a  series  of  maxima,  and  minima  at  constant  particle 
size,  their  number  being  .larger  the  wider  the  spectral  i-ange  used.  As 
shown  more  recently4  these  maxima  and  minima  move  to  longer  wavelengths 
and  become  shallower  as  the  syntem  becomes  more  and  more  heterodisperse . 
(See  particularly  Figure  5  end  6  in  ref.  4). 

In  order  to  correlate  the  spectral  location  of  6 -maxima  and 
minima  and  the  amplitude  of  the  spectral  oscillations  of  <f,  with  size 
distribution  curves,  a  certain  basic  type  of  distribution  curve  must  be 
assumed  unless  the  accessible  spectral  range  1b  extraordinarily  wide  or 
unless  a  series  of  t additional  light  scattering  criteria  are  used.  In¬ 
stead  of  using  any  of  the-  various  possible  well-known  distributions, 
such  as  Gaussian,  log  normal,  Maxwellian,  a  special  type  of  distribution 
woe  assumed  in  the  present  work:  a  distribution  which  fitted  veil  those 
found  in  emulsions.  Tie  reason  for  this  is  that  colloidal  emulsions 
roost  likely  will  rcpxesent  one  of  the  most  important  areas  where  this 
new  technique  is  va  riable  since  electron  microscopy  cannot  be  applied 
here  except  under  vezy  special  circumstances .  The  distribution  function 
selected  is 

?(r)  -  (r-rQ)  exp  £-  f(r-rQ)  /s]3J,  r^rQ,  (1) 

*  °»  r C^- 

where  Cf(r)  dr  is  the  number  of  particles  per  unit  volume  vith  radii 
between  r  and  r  +  dr,  C  So  a  normalisation  constant,  and  rQ  and  s  are 
parameters  cheractoxlzing  the  distribution;  r  is  the  radius  of  the 


6.  <A~  2 whsxo^s a  the  radius  of  the  sphere  and^io  the  v«ave« 

length  in  the  surrounding  medium. 

7.  W.  Heller,  W.  -J.  Pangonio,  and  N.  A.  Econonou,  J.  Chen.  Fhys. 

971  (1961). 
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smllast  particle  present  in  consequential  numbers  while  a  determines 

the  nodal  radius  r  ,  the  half  vidth  v,  and  the  "half  spread"  (r -r  ) 
n  no' 

of  the  distribution  through  the  relations 


w  -  O.9015  3 


r  -r 

31  0 


,-l/2o 
0  ‘  8 


(2) 

(3) 


The  type  of  distribution  curve  represented  by  eq.  (l)  1b  similar  to  a 
log-normal  distribution  inasmuch  as  it  is  positively  shoved  as  generally 
found  in  colloidal  oyetemo.  It  differs  from  it  and  from  any  other  dis¬ 
tribution  curve  by  the  fact  that  it  generally  begins  at  a  finite  <£  (at 
a  finite  particle  site).  Shis  has  two  reasons:  First,  particle  else 
distributions  within  the  Rayleigh  range  cannot  be  determined 

without  additional  criteria,  the  respective  section  of  an  optical  dis¬ 
tribution  curve  being  therefore  meaningless  irrespective  of  the  type  of 
distribution  assumed.  Secondly,  and  most  Importantly,  the  smallest 
limiting  particle  size  determined  with  this  type  of  function  represents 
well  the  number  of  particles  present  in  consequential  numbers.  This 
therefor©  introduces  an  additional  practically  useful  parameter  for  the 
characterization  of  perticle  size  distributions. 

I21  hetorodisperae  systems,  the  scattering  ratio  is 

6*  '*J* J  //  Cf(r',  Cf(r)dr  (4) 

Ears  J  is  intensity  of  light  scattered,  at  an  angle  of  90°  with  respect 

to  the  incident  beam,  by  a  single  particle  of  radius  r,  per  unit  solid 
angle  and  unit  intensity  of  the  incident  beam.  On  changing  the  variable 
from  r  to  <*. 

C'-  F//(p,q)  /F^(p,q).  (5) 

The  quantities  p  and  q  are  defined  by 

P  -  2r  tt(X  ,  (6) 


4  -  2s  vi  X' 


(7) 
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The  <T -values  corresponding  to  all  p  and  q  values  that  say  practically 
be  encountered  in  colloidal  dispersions  of  non-absorbing  spheres  have 

ft 

been  tabulated  for  the  relative  refractive  indices  m  -  1.05  (0-05)  1.J0 
using  ao  a  theoretical  basis  tho  Mie -theory^ . 

Ibr  practical  work  it  is  convenient  to  normalize  the  data  with 
respect  to  the  green  mercury  line  as  the  reference  vavo length,  In 

the  present  work,  where  the  dispersing  medium  is  water  at  25°^  ^  =  4093 . 57A  • 


The  corresponding  normalized  p  and  q  values 

»R  ■  a  TT*J  l  B  (8) 

”-8’  28  *IX*  <9) 

d 

As  already  shown7  the  practically  inportant  quantities 

w  =  0.05873q-j,  (30) 

Vr0  =  °*04517qR,  (11) 

rn  «  O.Ob517^  +  0.06515PR,  (12) 

ro  «  0.06515pr  (i*3) 


Tbs  actual  procedure  in  deriving  distribution  curves  consists  of 
the  following  steps:  (l)  spectra  of  the  scattering  ratio  are  obtained  at 
various  concentrations;  (2)  the  spectrum  pertinent  to  infinite  dilution  is 
derived;  (3)  a  few  exploratory  theoretical  spectra  are  derived  for  several 
pairs  of  pR  and  ^-values;  (4)  the  exact  pR  and  ^  values  are  decided  upon 


8.'*^i^»i»here^2_  is  the  refractive  index  of  the  spheres  and ✓4'.  that  of 
the  medium.  * 

9-  Ttebles  of  Scattering  Functions  for  Heterodlroerae  Systems. 

A.  P,  Stevenson  and  W.  Heller;  flaync  State  University  Press;  Detroit, 
Michigan,  1961. 
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which  vj II  give  the  best  possible  fit  of  the  theoretical  and  experi¬ 
mental  anectrum.  A  number  of  variations  of  this  basic  principle  are  poa- 

g 

oible  and  have  been  discussedT. 

In  tbe  present  work,  the  theoretical  spectra  giving  the  best  fit 
■•('re  obtained  by  means  of  hand  calculations,  in  order  to  get  as  intimate 
a  knowledge  of  the  sensitivity  of  the  theoretical  spectra  to  p^  and 
aa  possible.  In  systematic  practical  application  of  the  method,  the  fit¬ 
ting  procedure  by  means  of  electronic  computer  is  obviously  indicated  on 
account  of  the  considerable  amount  of  time  that  thus  can  be  saved.  Steps 
(3)  and  are  then  reduced  to  a  single  operation. 

III.  EXPEEIMEilTAL 

1.  The  Syr  teas  Investigated 

Ihc  experimental  tests  were  carried  out  on  heterodisperse 
latices  of  polystyrene,  obtained  by  selectively  mixing  portions  of  eighteen 
relatively  monc-disperse  latices.  The  mean  particle  size  of  the  monodisperoe 
Inti coo  varied  cyctematicalJy  from  753  to  1,300  m^y10  in  intervals  of  about 
30  m,.y.  Therefore  smooth  distribution  curves  of  almost  any  type  could  be 
obtained  by  proper  nixing.  The  polystyrene  systems-in  which  m  closely 
approximates  l.J-X)  at  5^1A-vere  attractive  also  because  the  light  scatter¬ 
ing  of  nonodi sperre  polystyrene  latices  had  already  been  investigated 
extensively^'.  Three  beterodieperse  systems  were  thus  prepared,  one  in 
which  the  distribution  approximated  that  given  by  Equation  1  (H.D.I.),  a 
second  ir  which  the  distribution  was  antioynmetrical  i.e.  negatively 
skewed,  (H.D.  2.)  and  a  third  which  was  approximately  Gaussian,  (H.D.  3)* 

Each  of  the  3  prepared  stock  mixtures  had  a  concentration  of  10  g 
of  polymer  per  1,000  g  latex.  Prom  each  of  them  a  series  of  dilution 
were  made  so  that  tbe  -spectra  could  bo  extrapolated  to  zero  concen- 


10.  The  monodisperoe  samples  vere  kindly  provided  by  Dr.  J.  W.  Vanderhoff 
of  the  Dow  Chemical  Company. 

11.  R.  M.  Ittbibian,  W.  Heller,  and  J.  H.  Epel,  J.  Coll.  Sci.  U,  195  (1956) 
and  later  papers,  ibid. 
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t  rat  ion.  Hvs  highest  and  lowest  concentration  optically  investi¬ 
gated  within  these  three  dilution  aeries  were  <^1.4o  x  10  and 
(5^8)  x  lo'3  g  of  polymer  per  1,000  g  of  latex.  For  technical  details, 
reference  may  be  made  to  the  Investigations  on  mono&i aperae  systems'3 
as  nee  they  were  identical  except  for  the  following:  Centrifuging  was 
omitted  and,  prior  to  filtering,  stabilizer  was  added  and  ultrasonic 
irradiation  carried  out  for  10  minutes..  This  was  found  to  effectively 
disperse  any  aggregates  that  may  have  been  present. 


2 ,  Ifce  Optical  Measurements  , 

The  apparatus  described  previously3  was  used  for  the 
present  work  also  after  applying  the  following  modifications:  (1)  A  150 
candle  power  a.c,  Pbintollt.e  arc  lamp  was  employed  as  an  incandescent 
light  source  since  it  combined  high  brightness *at  all  wavelengths  in  the 
visible  with  virtually  constant  outp-5t  (2)  Among  the  various  diaphragms 
used  in  the  apparatus  (See  Fig.  1  In  reference  3)  the  following  were 
increased  in  aperture:  A-  (4-*  5-4  mm.),  Dg  (U.O-^5.2  am.}, 

D...  (4.1  •♦6.4  nan.).  As  done  previously  (6.7  nut. )  was  generally 
omitted .  "*  ’ 


Tha  spectra  wore  investigated  between  the  limiting  wavelengths 
of  A510  and  5970A  for  H.D.l  and  2*.  and  between  4540  and  6000A  for 
B.D.3.  In  order  to  increase  the  intensity  of  the  incident  beam,  a 
slightly  larger  entrpace  slit  to  the  monochromator  A  (see  fig.  I3)  was 
used  than  previously^.  The  resulting  increase  in  the  low  intensity 
polychromatic  background  of  the  narrow  spectral  band  emerging  from  the 
monochromator  was  largely  compensated,  however,  by  passing  the  light 
bean,  prior  to  its^entry  into  the  monochromator,  through  colored 
aqueous  solutions  .  The  width  of  the  spectral  band  obtained  from  the 
Qonechroraatcr  and.  its  variation  with  wavelength  will  be  taken  into  ac¬ 
count  in  the  results  tc  be  discussed  be low A- . 


12.  H»bc  solutions,  acidified  if  necessary,  contained  two  or  more 
ceuponents  of  the  following  salts  in  a  suitable  concentration: 
Ctt'HD^,  CTN03)3.  CuCl?,  KgCr^,  Co(3C^),  and  lttClg. 

13  The  width  of  the  spectral  band  entering  the  apparatus  ■varied  be¬ 
tween  the  blue  and  the  orange  from  11  to  32,  5  to  12,  and  6 
to  16  millimicrons  in  the  experiments  with  H.D.l,  H.D.2,  and  H.D.3, 
respectively.  The  more  important  half  width  of  the  intensity  dis¬ 
tribution  throughout  the  spectral  band,  which  is  decisive  for  the 
response  of  the  photocell,  was  of  course  considerably  smaller,  but 
it  is  the  former  which  will  be  used  in  determining  the  m»ytt»nn 
uncertainty  of  the  results. 


The  solid  an^le  •*  controlled  by  the  position  of  diaphragm  DIP 
was  varied  within  a  relatively  narrow  range  depending  on  the  scattered 
intensity,  required  for  a  measurement .  Tt  had  a  minimum  value  of 
1.23x  10  3  and  a  maximum  value  ot  9.6  x  10~3  ete radians .  The  latter 
was  used  in  the  majority  of  measurements .  This  compares  with  the 
standard  solid  angle  of  1.5  x  10 “3  steradians  used  previously.  This 
angle  was  still  within-, the  range  where  the  value  of  the  scattering  ratio 
is  not  affected  as  yet3. 

The  anisotropy  of  the  photomultiplier  tube  was  determined 
quantitatively  by  Bpecia^  experiments  and  was  taken  into  account  in 
the  numerical  data  of  <f1  ■ 

The  actual  determinations  of  the  scattering  ratio  were  carried 
out  according  to  the  detailed  technique  described  previously3. 

3  Electron  Microscopy 

In  view  of  the  mixing  procedure  used,  the  synthetic  sire 
distribution  curve  was,  of  course,  a  priori  known,  nevertheless,  an 
electron  microscopic  analysis  of  the  final  mixtures  was  also  carried 
out.  The  optical^atc-  to  be  reported  below  are  checked  against  these 
alternate  results  The  analysis  of  the  electron  photomicrographs  was 
carried  out  on  their  projections  onto  a  large  screen  giving  an  overall 
magnification  of  1^0,000 .  Only  the  central  portions  of  the  individual 
micrographs  were  considered  in  order  to  eliminate  the  effect  of  minor 
distortions  in  the  peripheral  sections  which  may  not  necessarily  be 
apparent  visually.  In  order  to  be  absolutely  certain  of  the  results 


lb.  The  anisotropy  may  be  expressed  in  terms  of  d  ///da  .  Here  d 
is  the  actual  deflection  of  the  galvanometer  in  the  absence  of 
any  scattering  solution  with  the  photomultiplier  tube  mounted  with 
its  long  axis  in  the  plane  of  observation,  its  customary  position 
in  the  apparatus  used.  The  subscripts  //  anda  pertain  to  an 
incident  polarized  beam  vibrating  parallel  and  perpendicular, 
respectively,  to  the  plane  of  observation.  The  ratio  varied  fr cm 
0.96b  to  0.971  between  4500  and  6000  A.  viz.,  by  only  O.TJl  over 
the  entire  spectrum  used. 

15-  The  electron  photomicrographs  were  kindly  provided  by  Dr. 

J.  B.  L.  Watson  of  the  Edsel  B.  Ford  Institute  for  Medical 
Research,  Detroit,  Michigan . 
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for  H.D.l  and  H.D.2,  all  particles  in  the  central  section  considered 
were,  in  addition,  measured  in  two  orthogonal  directions  and  all  those  - 
results  were  rejected  in  which  the  axial  ratio  differed  from  unity  by 
more  than  5$- 

IV.  RESULTS 

1.  Distribution  Curve  of  the  Type  Given  by  Bqa.  1. 

The  small  black  discs  in  Fig.  1  represent  the  most 
probable  values  of  the  cf'-data  obtained  experimentally  on  H.D.l.  The 
rectangle  within  which  the  most  probable  value  ic  contained  or  horizontal 
line  which  passes  through  it  give  the  extreme  limits  in  the  uncertainty 
of  these  date.  The  exact  width  of  the  band  of  incident  radiation^  is 
indicated  by  the  horizontal  dimension  of  the  rectangle  (length  of  the 
line)  while  the  maximal  uncertainties  in  the  numerical  value  of 
(I  //  j’Lj_  )o  (The  subscript  o  signifies  extrapolation  to  zero  concen¬ 
tration)  are  represented  by  the  height  of  the  rectangle.  Among  all  pos¬ 
sible  theoretical  curves  which  fit  the  experimental  data  on  varying  p^ 
and  in  steps  of  0.2,  the  choice  could  be  narrowed  down  to  the  set  of  four 
lightly  dravn  or  dashed  curves  given.  These  curves  were  calculated  from 
Eq.  6  in  Ref.  jji^ .  They  pertain  to  the  parameters  (p^,  q^)  of  (6. 8, 1.4), 
(6.8,1 .6),  (7.0, 1.4),  (7.0,1. 6).  interpolation  between  the  four  curves 
one  can  find  one  which  fits  the  experimental  data  most  closely;  this  curve, 
corresponding  to  pR  ■  6.84,  q^  «  1.49,  is  represented  by  Curve  I  in  Fig.  1. 
It  gives  the  fully  drawn  distribution  curve  in  Fig.  2.  (The  ordinate  Cf(r) 
is,  as  stated,  proportional  to  the  number  of  particles,  per  unit  volume, 
whose  radius  is  r.) 


16.  The  theoretical  data  were  calculated  assuming  that  n  »  1.20  over 
the  entire  spectral  range.  The  dispersion  of  m,  which  be 

taken  into  account  when  significant,  was  found  not  to  alter  the 
general  pattern  of  the  results  here  and  could  therefore  be 
neglected. 
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The  electron  microscopically  obtained  histogram  of  the  distri¬ 
bution  curve  ^analysis  of  982  particles) is  also  shown  in  Fig.  217.  Xu 
view  of  the  expanded  scale  of  the  abscissa;  the  agreement  must  be  con¬ 
sidered  as  very  satisfactory.  The  peal;  of  the  distribution  curve 
representing  the  uodal  diameter  (Dffl)  as  obtained  by  the  two  methods 


differs  by  5Jk,  the  half  width,  by  2.2$,  the  "half  spread"  2(r  -r  )  by  only 

SI  o 

Z$  and  the  aaallest  particle  diameter  DQ  by  6$.  Tbs  theoretical  curve  In 


Fig.  1  which  satisfied  the  experimental  data  the  least  (pR  «  7.0,  q^  »  1.6) 

gives  tbs  dashed  distribution  curve  XX  in  Fig.  2.  It  is  readily  seen  that 

it  does  not  differ  ouch  from  the  distribution  curve  obtained  by  inter- 

ifl 

point  ion.  fiance,  one  may  dispense  with  Interpolation  if  optimum  accu¬ 
racy  is  not  essential. 


2.  Negatively  Skewed  Distribution  Curve 

The  experimental  data  tor  H .  D.  2^4iatribution  curve  aoti- 
symnetrical  to  that  postulated  by  Sq.  l)are  given  in  Fig.  3-  Die  theo¬ 
retical  curw  which  cooes  closest  to  an  acceptable  fit  with  the  exper¬ 
imental  data  le  obtained  by  interpolation  between  4  pairs  of  pR  and  q^ 
values.  H>e  parameters  of  this  Curve  (I),  pR  «*  8.15,  ^  -  O.85,  yield 
the  normalised  distribution  curve  1  in  Fig.  4.  Cccparison  with  the 
electron  nicroscoplc  histogram  shows  that  the  number  distribution  ob¬ 
tained  is  quantitatively  not  satisfactory.  It  Is  noteworthy,  however, 
that  the  derived  distribution  function  does,  even  in  case,  repro¬ 
duce  the  electron  microscopic  modal  diameter  to  within  4$. 


17 .  Doc  optical  distribution  curve  is  normalized  with  respect  to  the 

election  microscopic  histogram  00  that  they  each  subtend  equal  areas. 

•  This  procedure  was  adopted  with  all  the  systems  Investigated. 

l£.  Normalization  of  both  experimental  end  theoretical  data  repre¬ 
sented  in  Fig.  1  with  respect  to  tbs  data  at  tbs  reference  wave¬ 
length  may  facilitate  the  analysis  and  in  favorable  cases  improve 
the  accuracy  of  the  results. 


Xus  fact  that  the  distribution  curve  obtained  from  light 
scattering  deviated  bore  from  the  actual  one  is  of  courao  not  sur¬ 
prising.  Indeed,  it  is  very  gratifying  to  see  from  Fig*  3  ae  compared 
to  Fig.  1  that  a  radical  deviation  of  the  distribution  curve  from  that 
assumed  cxmiffcuts  ltoelf  immediately  by  the  fact  that  no  theoretical 
suectrun  can  be  ftound  which  satisfies  the  experimental  one  in  satis¬ 
factory  approacinptico  over  a  sufficiently  extensive  spoctral  range. 

She  risk  of  obtaining  fictitious  distribution  curves  from  light  scat¬ 
tering  spectra  is  therefore  quite  small  since  the  impossibility  of  a 
satisfactory  fit  between  theoretical  and  experimented  spectra  represents 
an  automatic  check. 

The  question  naturally  arises  as  to  what  may  usefully  be  done  in 
absence  of  a  satisfactory  fit.  The  possibilities  of  resolving  the  pro¬ 
blem  will  be  discussed  in  section  V  with  emphasis  upon  one  possibility 
adopted  in  the  present  work. 

3 .  Gaussian  Distribution 

Considering,  finally,  H.D.3  with  its  normal  distribution 
curve,  Jig.  5  gives  the  experimental  of  )  data  and,  in  addition,  the 
interpolated  beet  fitting  theoretical /*(^  ) -curve  characterized  by 
pp  »  5.63,  qg  *  I.90.  The  size  distribution  obtained  from  this  best 
f  )  "Curve  is  compared  with  the  electron  microscopic  histogram  in 
Fig.  6.  The  modal  diameters  agree  to  within  5$  which  is,  like  in  the 
case  of  H.D.l  and  2,  within  tbs  uncertainty  of  the  electron  microscopic 
histogram.  The  half  width  agrees  to  within  27$,  and  the  value  for  the 
emallti&t  particles  present  in  consequential  numbers  to  within  l!t£. 

Normal  distribution  curves  can,  therefore,  be  analysed  satisfactorily 
with  the  method  outlined.  The  deviations  of  the  best  fitting  theoretical 
spectaxm  from  the  experimental  spectrum  are  distinct,  but  considerably 
smaller  than  in  the  preceding  instance  (H.D.2).  Therefore,  very  precise 
experimental  data  are  nooeseary  in  order  to  give  a  clear  lndicatlon-as 
apparent  from  Fig.  5 -that  the  distribution  curve  is  not  fully  in  line 
with  that  assumed. 


V.  DETERMINATION  OP  SIZE  DISTRIBUTION  CURVES  VEEN  THEORETICAL 
AND  IXPERIMEIfEAL  P'- SPECTRA  DO  NOT  MATCH  SATISFACTORILY . 


In  absence  of  a  satisfactory  fit  between  theoretical 

and  experimental  S -spectra -which  indicates  non-compliance  of  the  actual 

distribution  curve  with  that  assumed,  one  has  the  choice  between  four 

k 

possibilities,  three  of  which  have  already  been  briefly  discussed  . 

(l)  A  first  approximation  method  based  on  graphical  interpolation. 


The  sise  distribution  curves  for  H.D.l,  H.D.2,  and  H.D.3  were 
obtained  by  interpolation  for  the  proper  and  values  from  four  of 
those  nX)  curves  which  come  closest  to  satisfying  the  ejqperimental 


data.  An  alternate  method  since  then  adopted  consists  of  picking  two 
curves  only,  selected  in  such  a  manner  that  all  experimental  points  are, 
throughout  the  entire  spectrum,  contained  within  the  area  between  the 
two  curves.  One  thus  establishes  two  distribution  curves.  While  they 


differ  very  little  if  the  actual  distribution  curve  is  of  the  positively 


skewed  type,  i.e.  if  the  area  between  the  two  theoretical  spectra  la 
very  small,  an  appreciable  difference  is  found  In  other  circumstances. 


Graphical  interpolation  between  the  two  curves  yields  then  a  distri¬ 
bution  curve  approximating  that  of  the  system.  This  interpolated  curve 
clearly  deviates  from  the  type  given  by  eq.  (1).  It  may,  in  extreme  cases 
be  primarily  negatively  skewed  or  even  bimodal^. 


(2)  /.-second  approximation  method  based  on  the  use  of  a  two  term  equation. 

A a  already  pointed  out  previously^,  the  rigidity  with  which  the 
basic  type  of  distribution  is  fixed  by  eq.  (l)  can  be  removed  by  a 


19.  Since  in  practical  large  scale  application  of  the  methods  outlined 
hare,  computer  fitting  of  theoretical  and  experimental  spectra  is 
highly  desirable;  the  procedure  to  be  adopted  In  the  present  instance 
la  the  independent  match  of  two  halfs  of  the  spectrum  and  a  subsequent 
match  of  the  entire  spectrum.  If  the  minimal  sum  of  the  absolute 
deviations  (or  of  their  squares)  from  the  experimental  data  has  in  all 
three  Instances  approximately  the  same  value,  eq.  (1)  applies.  Other¬ 
wise,  the  computer  will  have  to  he  ordered  to  derive  from  the  two 
complete  spectral  curves  the  arithmetic  mean.  The  latter  result  is 
that  corresponding  to  the  first  approximation  method  outlined. 
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second  term.  'Phis  gives 


g(r)  *  f<r,ro,e)  +•  (C/C^)  f(r,r^,s‘)  <» 

*  (r-rQ)  exp  f>\  [r-ro  J +  (Cg/^)  'r-r^)  exp  (-[{r-i^)/a']^ 


Ibis  equation,  was  applied  to  H.D.2,  using  analytical  principles 
4 

discussed  previously  .  The  heavily  drawn  curve  II  in  Fig.  3  shows  the 
theoretical  <f  { X  )  curve  thus  obtained  and  the  heavily  drawn  curve  H  In 
Figo  4  represents  the  resulting  distribution,  It  agrees  satisfactorily 
with  the  electron  microscopic  histogram.  In  particular,  the  number  dis¬ 
tribution  is  essentially  correct  and,  in  addition,  the  tail  towards  the 
small  particle  size  range  is  reproduced  satisfactorily.  It  will  be  noted 
that  the  nodal  diereter  is  the  some  as  before  ,ron  using  rather  eq.  (l)]« 

Ibe  only  unavoidable  artifice  introduced  by  the  second  tern  is  the  sudden 
change  in  slope  whrch  occurs  in  the  pjescnt  instance  between  Dq  and  Dffi. 

In  practical  application,  it  would  be  reasonable  to  3000th  out  such  a  kink. 

On  using  the  two  term  equation,  it  is  necessary  to  know  the 
individual  theoret-.cal  values  of  1  //  and  I j  or  KI  //  and  KI  where  K  is  a 
constant  at  a  Riven  vavelength.  This  Information  is  included  in  the  tables 

9 

referred  to  above  , 

The  task  of  finding  the  proper  and  values  can  be  simplified 
considerably  by  using  the  two  limiting  distribution  curves  established  ac¬ 
cording  to  the  first  approximation  procedure  discussed  above  (V,  l) .  Any 
of  the  two  curves  provides  the  Vp  and  q^  values  while  the  other  provides 
the  and  q^  values.  The  proper  C^/C^  is  found  by  successive  approximation, 
but,  in  general,  no  error  of  consequence  will,  be  committed  on  assuming  that 
Cg/C^  =  1.0. 

(3)  Other  Possibilities 

Two  further  possibilities,  which  already  have  been  discussed 
briefly  are  at  present  under  investigation,  but  no  statement  can  be  made 
as  yet  on  their  practical  advantage  over  the  2  preceding  procedures. 

first  of  all,  one  may,  in  absence  of  a  satisfactory  fit,  test  the  ex- 
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perlzrental  (/-spectrum  successively  against  tbsoretioal  spectra  result* 
lag  from  2  or  3  characteristically  different  theoretical  prototypes  of 
dlotritoutlon  fUnctlouo.  Silo  procedure  makes  the  use  of  an  electronic 
computer  mandatory. 

The  second  additional  possibility  whlch-tbeoretloally  at  least* 
is  the  most  attractive  one  consists  of  deriving  from  scattering  data  a 
distribution  curve  without  making  any  assuqption  whatsoever  about  its 
character.  While  this  is,  In  principle  possible,  It  is  nsceasazy  to 
increase  the  number  of  experimental  data  considerably  if  a  single  val¬ 
ued  answer  is  to  be  obtained.  In  other  words,  in  addition  to -spectra, 
turbidity  spectra  (discussed  in  a  subsequent  paper)  and  the  variation  of 
scattering  with  the  angle  of  observation  Bust  be  included  as  erqperlmental 
arguments. 

VI.  Simplification  of  Procedure  If  the  Objective  is  Merely 
The  Determination  of  the  Peak  of  The  Distribution  Curve 
Or  of  the  Member  Average  Particle  Size. 

Fbr  many  purposes,  particularly  in  industrial  control 
work,  it  is  less  important  to  know  all  the  details  of  a  particle  size 
distribution  than  it  la  to  determine  the  particle  diameter,  D  •  83?  at 

0  IB 

which  the  distribution  has  its  peak  (modal  diameter)  or  the  number 
average  particle  disaster,  Dq.  The  simplifications  which  are  then  pos¬ 
sible  In  the  procednrer.  described  become  apparent  on  the  basis  of  jMMLa  T- 
It  lists  in  the  first  row  Dq  (=2rQ)  and  as  obtained  tor  all  three 
systems  from  (/-spectra  using  eqn.  (l) .  The  percent  deviation  of  tbs  re¬ 
sults  relative  to  electron  microscopic  determination  is  given  in  paren¬ 
theses.  (The  deviation  is  considered  positive  if  the  light  Boattaring 
value  is  larger)  .  The  second  row  oontains  the  corresponding  data  ob¬ 
tained  on  using  tbs  two  term  .-go.  (lb-)  tor  H.D.2. 

It  is  seen  that  the  deviation  in  is  throughout  within  the  limits 
of  uncertainty  of  electron  microscopic  measurements  (about  5£)  even  on 


20.  The  deviations  were  evaluated  by  means  of  polygons  inscribed  in  the 
histograms  of  Jig's.  2,  k,  and  6.  They  were  exaead  prior  to  re¬ 
production  of  these  figures. 


o 


applying  eqn,  (l)  to  a  system  with  a  characteristically  different  dis¬ 
tribution  curve  (B.D.2).  A  very  rough  fit  of  ejqperimantal  and  theoret¬ 
ical  <f-epectra  la  therefore  sufficient  in  order  to  obtain  a  satis¬ 
factory  value  fbr  D&.  This  Is  airily  due  to  the  ftxct  that 

-  0.09303^  ^  +  0.13030  pjj  (15) 

i.e.  Dffl  varies  only  slowly  with  and  pR.  Sun,  the  four  distribution 
curves  derived  from  the  four  theoretical  spectra  in  Pig.  1  differ  in 
by  less  than  2f>.  Pile  fact  and,  particularly,  the  very  satisfactory 
Dm  value  Obtained  for  H.D.2  shorn  also  clearly  that  coincidence  of  DQ 
obtained  both  by  light  scattering  and  by  another  method  is  no  guarantee 
whatsoever  that  the  distribution  curve  assumed  applies  to  the  system 
under  investigation.  It  is  important  to  note  that  this  applies,  of 
course,  also  to  other  possible  experimental  criteria,  such  as  turbidity 
spectra  or  angular  variation  of  scattering  and  to  the  use  of  other 
theoretical  distribution  functions.  In  contradistinction  to  1^,  the 
value  of  Dq  is  very  sensitive  to  the  type  of  distribution  assumed  as  ap¬ 
parent  fnaa  the  data  of  B.D.2  and  B.D.3  as  compared  to  H.D.l21. 

The  statements  mode  with  regard  to  Dn  apply,  in  a  first  approxi¬ 
mation,  also  to  the  nuatoor  avenge  particle  diameter, 

5n  -  0.06515  (1.476  ^  ♦  8.000  jjj)  (16) 

as  the  siisUority  of  tbo  numerical  coefficients  in  eqpa.  (15)  and  (l6) 
indicates. 

A  particularly  staple  situation  arises  If  is  appreciably  larger 
than  qa.  It  follows  then  from  e«i.  (15)  sod  (16)  that,  in  a  first 
•Bproxlmstlon  DB  and  DQ  respectively  axe  directly  proportional  to  p^. 

In  other  words,  it  is  then  possible  to  derive  Da  (5q)  from  a  slt«le 
aesmirwwat  at  the  reference  wavelength.  Tt&»  Is  tasted  in  the  third 
rear  of  Table  I  for  D^.  These  data  are  sisply  Obtained  from  measurements 


21.  3ha  electron  microscopic  value  of  r  was  obtained  by  disoouatlam 
those  ami  last  particles  whose  total  amber  was  loss  tbn  1  per  cent 
of  all  the  particles.  . 
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at  the  green  Bg  -  line  on  associating  with  the  t? -value  that  diameter 
which  would  obtain  if  the  system  were  onsodlsperse .  (Sbr  this  purpose 
the  theoretical  data  given  In  reference  7  were  used).  The  error 
relative  to  the  rigorous  D  -values  in  row  one  of  Table  I  is  clearly 
4-5 throughout.  Consequently,  this  procedure  is  quantitatively  acceptable 
up  to  a  q^/p^  ratio  of  about  0.35 •  This  means,  in  practical  terms,  a 
distribution  curve  whose  half  width  is  not  more  than  about  one  third 
as  large  as  the  smallest  particles  present  in  consequential  numbers. 

Thus,  If,  for  instance,  a  distribution  starts  effectively  at  about 
0.5  microns,  Dm  (and  Dq)  will  be  obtained  from  a  single  monochromatic 
•measurement  with  an  error  not  in  excess  of  5$  if  the  haliVldths  is  not 
in  excess  of  160  millimicrons  *  This  is  essentially  equivalent  to  re¬ 
quiring  that  no  particles  should  be  present  in  consequential  numbers 
below  500  millimicrons  nor  above  800-8 50  millimicrons. 

If  is  appreciably  larger  than  p^,  and  Dq  will,  in  a  first 
approximation  be  proportional  to  q^,  again  according  to  eqs.  (15)  sad 
(16).  Therefore  in  strongly  heterodisperse  systems  in  which  signifi¬ 
cant  numbers  of  particles  are  present  even  within  the  lowest  section 
of  the  colloidal  range,  the  approximate  value  of  (5Q)  obtained  at 
a  single  wavelength  is,  in  first  approximation  a  relative  measure  of 
the  degree  of  heterodispersion  of  the  system. 

VH.  THE  BASSE  0?  PRACTICAL  USEFULHES8  OP  THE  I, TORT  SCATT3HBG  METHOD 
m  COSBECriOfl  WISH  THE  problem  op  size  distributiohs. 

One  of  the  most  attractive  features  of  size  distribution  deter¬ 
minations  from  light  scattering  is  the  fact  that  the  systems  under 
investigation  are  not  interfered  with,  in  contradistinction  to  all 
other  known  methods.  Prom  the  theoretical^  and  experimental,  investi¬ 
gation-including  experimental  evidence  obtained  after  conclusion  of 
oo 

the  present  work  it  follows  that  the  singlevaluedness  and  accuracy 
of  data  obtained  from  scattering  depends  on  the  spread  of  the  distri- 


22.  Data  obtained  in  this  laboratory  by  Rodney  Wu,  and,  more  recently, 
by  Bettye  Greene  (to  be  published  later). 


-  16  - 


button  curve  and  on  the  size  of  the  smallest  particles  present  in 
consequential  numbers.  The  role  of  the  latter  ia  relatively  minor. 

Prom  Fig.  2  (ref.  4),  it  follows  that  the  sensitivity  of  the  method 
should  increase  with  p.  On  using  visible  light,  an  optimum  can,  how¬ 
ever,  be  expected  et  a  particle  diameter  in  excess  of  several  microns 
since  then  the  number  of  maxima  and  minima  per  unit  spectral  range  be¬ 
comes  too  large  to  he  measured  conveniently.  The  width  of  the  distribution 
curve,  on  the  other  hand,  is  far  more/  important; Fig .  3  (ref.  5)  shows 
that  at  q-Yalxies*7''.5,^/  varies,  at  constant  p,  only  slowly  and  nearly 
linearly  with  q.  The  sensitivity  of  the  method  is  therefore  optimal 
at  q-values  <T2.5»  This  corresponds  to  a  distribution  curve  with  a  half 
width  not  in  excess  of  about  200  millimicrons,  and  an  entire  spread  not 

in  excess  of  about  400-900  millimicrons,  if  visible  light  is  used^. 

4 

It  follows  from  Fig.  3  that  for  at  a  given 

Aq  +  Bp  -  Cpq  ‘S' constant  (17) 

so  that  in  the  absence  of  a  numerically  large  spectral  variation  of  {T , 
(maxima  and  minima  being  then  essentially  absent),  a  relatively  large 
number  of  p  and  q  combinations  may  satisfy  the  same  spectrum.  This 
problem  Is,  clearly,  not  limited,  to  the  use  of  eq.  (l)  nor  to  the  use 
of  the  experimental  criterion  used  here  (<f -spectra).  It  is  bound  to 
be  encountered  also  on  basing  the  analysis  on  other  types  of  distri¬ 
bution  or  on  using  aa  experimental  criteria  turbidity  spectra  or  the 
variation  of  scattering  with  the  angle  of  observation.  The  lack  of 
singlevaluedness  of  results  to  be  expected,  therefore,  in  very  hetero- 
dlsperse  systems  can  be  relieved,  however,  by  using  in  addition  to 
-spectra  the  two  essentially  equivalent  criteria  Just  named.  A 
subseqient  paper  in  this  series  will  be  concerned  with  this  possibility. 


23*  The  allowed  spread  increases  in  direct  proportion  to  the  wave¬ 
length.  It  is  at  leant  1000  times  larger  at  radar  wavelengths. 
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Figure  1.  /j~- spectrum  of  H.D.l,  a  system  conforming  to  the 
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black  circles:  most  probable  experimental  value, 
rectangles:  maximum  experimental  uncertainty, 
curves:  theoretical  spectra. 

Figure  2.  Size  distribution  curves  of  H.D.l. 

I.  Derived  from  best  fitting  interpolated  theoretical 
& -spectrum. 

II.  Derived  from  approximately  fitting  theoretical 
d'-spectrua. 

Histogram:  electron  microscopic  distribution  curve. 

Figure  3.  ^-spectrum  of  H.D.g  (negatively  skewed  distribution). 

I.  Best  possible  theoretical  cT-spectrun  achieved 
by  means  of  eq.  (l) . 

II.  Better  theoretical,  f  -spectrum  obtained  on  using 
equation  (l4). 

h.  Size  distribution  curve  of  H.D.2. 


I.  Using  one-term  equation  (l). 
XI.  Using  two-term  equation  (14). 


Figure  5. 


(Gaussian  distribution) . 


Figure  6.  Size  distribution  of  H.D.3. 
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